Search for: All records

Abstract In this paper, we present an open-source software library that can be used to numerically simulate the advection and diffusion of a chemical concentration or heat density in a viscous fluid where a moving, elastic boundary drives the fluid and acts as a source or sink. The fully-coupled fluid-structure interaction problem of an elastic boundary in a viscous fluid is solved using Peskin’s immersed boundary method. The addition or removal of the concentration or heat density from the boundary is solved using an immersed boundary-like approach in which the concentration is spread from the immersed boundary to the fluid using a regularized delta function. The concentration or density over time is then described by the advection-diffusion equation and numerically solved. This functionality has been added to our software library, IB2d , which provides an easy-to-use immersed boundary method in two dimensions with full implementations in MATLAB and Python. We provide four examples that illustrate the usefulness of the method. A simple rubber band that resists stretching and absorbs and releases a chemical concentration is simulated as a first example. Complete convergence results are presented for this benchmark case. Three more biological examples are presented: (1) an oscillating row of cylinders, representative of an idealized appendage used for filter-feeding or sniffing, (2) an oscillating plate in a background flow is considered to study the case of heat dissipation in a vibrating leaf, and (3) a simplified model of a pulsing soft coral where carbon dioxide is taken up and oxygen is released as a byproduct from the moving tentacles. This method is applicable to a broad range of problems in the life sciences, including chemical sensing by antennae, heat dissipation in plants and other structures, the advection-diffusion of morphogens during development, filter-feeding by marine organisms, and the release of waste products from organisms in flows. 

Synopsis Computational models of aquatic locomotion range from modest individual simple swimmers in 2D to sophisticated 3D multi-swimmer models that attempt to parse collective behavioral dynamics. Each of these models contain a multitude of model input parameters to which its outputs are inherently dependent, that is, various performance metrics. In this work, the swimming performance’s sensitivity to parameters is investigated for an idealized, simple anguilliform swimming model in 2D. The swimmer considered here propagates forward by dynamically varying its body curvature, similar to motion of a Caenorhabditis elegans. The parameter sensitivities were explored with respect to the fluid scale (Reynolds number), stroke (undulation) frequency, as well as a kinematic parameter controlling the velocity and acceleration of each upstroke and downstroke. The input Reynolds number and stroke frequencies sampled were from [450, 2200] and [1, 3] Hz, respectively. In total, 5000 fluid–structure interaction simulations were performed, each with a unique parameter combination selected via a Sobol sequence, in order to conduct global sensitivity analysis. Results indicate that the swimmer’s performance is most sensitive to variations in its stroke frequency. Trends in swimming performance were discovered by projecting the performance data onto particular 2D subspaces. Pareto-like optimal fronts were identified. This work is a natural extension of the parameter explorations of the same model from Battista in 2020. 

Synopsis Computational scientists have investigated swimming performance across a multitude of different systems for decades. Most models depend on numerous model input parameters and performance is sensitive to those parameters. In this article, parameter subspaces are qualitatively identified in which there exists enhanced swimming performance for an idealized, simple swimming model that resembles a Caenorhabditis elegans, an organism that exhibits an anguilliform mode of locomotion. The computational model uses the immersed boundary method to solve the fluid-interaction system. The 1D swimmer propagates itself forward by dynamically changing its preferred body curvature. Observations indicate that the swimmer’s performance appears more sensitive to fluid scale and stroke frequency, rather than variations in the velocity and acceleration of either its upstroke or downstroke as a whole. Pareto-like optimal fronts were also identified within the data for the cost of transport and swimming speed. While this methodology allows one to locate robust parameter subspaces for desired performance in a straight-forward manner, it comes at the cost of simulating orders of magnitude more simulations than traditional fluid–structure interaction studies. 

While pendulums have been around for millennia and have even managed to swing their way into undergraduate curricula, they still offer a breadth of complex dynamics to which some has still yet to have been untapped. To probe into the dynamics, we developed a computational fluid dynamics (CFD) model of a pendulum using the open-source fluid-structure interaction (FSI) software, IB2d. Beyond analyzing the angular displacements, speeds, and forces attained from the FSI model alone, we compared its dynamics to the canonical damped pendulum ordinary differential equation (ODE) model that is familiar to students. We only observed qualitative agreement after a few oscillation cycles, suggesting that there is enhanced fluid drag during our setup’s initial swing, not captured by the ODE’s linearly-proportional-velocity damping term, which arises from the Stokes Drag Law. Moreover, we were also able to investigate what otherwise could not have been explored using the ODE model, that is, the fluid’s response to a swinging pendulum—the system’s underlying fluid dynamics. 

Computational Fluid Dynamics (CFD) models are being rapidly integrated into applications across all sciences and engineering. CFD harnesses the power of computers to solve the equations of fluid dynamics, which otherwise cannot be solved analytically except for very particular cases. Numerical solutions can be interpreted through traditional quantitative techniques as well as visually through qualitative snapshots of the flow data. As pictures are worth a thousand words, in many cases such visualizations are invaluable for understanding the fluid system. Unfortunately, vast mathematical knowledge is required to develop one’s own CFD software and commercial software options are expensive and thereby may be inaccessible to many potential practitioners. To that extent, CFD materials specifically designed for undergraduate education are limited. Here we provide an open-source repository, which contains numerous popular fluid solvers in 2 D (projection, spectral, and Lattice Boltzmann), with full implementations in both MATLAB and Python3. All output data is saved in the . v t k format, which can be visualized (and analyzed) with open-source visualization tools, such as VisIt or ParaView. Beyond the code, we also provide teaching resources, such as tutorials, flow snapshots, measurements, videos, and slides to streamline use of the software. 

Most biological functional systems are complex, and this complexity is a fundamental driver of diversity. Because input parameters interact in complex ways, a holistic understanding of functional systems is key to understanding how natural selection produces diversity. We present uncertainty quantification (UQ) as a quantitative analysis tool on computational models to study the interplay of complex systems and diversity. We investigate peristaltic pumping in a racetrack circulatory system using a computational model and analyse the impact of three input parameters (Womersley number, compression frequency, compression ratio) on flow and the energetic costs of circulation. We employed two models of peristalsis (one that allows elastic interactions between the heart tube and fluid and one that does not), to investigate the role of elastic interactions on model output. A computationally cheaper surrogate of the input parameter space was created with generalized polynomial chaos expansion to save computational resources. Sobol indices were then calculated based on the generalized polynomial chaos expansion and model output. We found that all flow metrics were highly sensitive to changes in compression ratio and insensitive to Womersley number and compression frequency, consistent across models of peristalsis. Elastic interactions changed the patterns of parameter sensitivity for energetic costs between the two models, revealing that elastic interactions are probably a key physical metric of peristalsis. The UQ analysis created two hypotheses regarding diversity: favouring high flow rates (where compression ratio is large and highly conserved) and minimizing energetic costs (which avoids combinations of high compression ratios, high frequencies and low Womersley numbers). 

Jellyfish are majestic, energy-efficient, and one of the oldest species that inhabit the oceans. It is perhaps the second item, their efficiency, that has captivated scientists for decades into investigating their locomotive behavior. Yet, no one has specifically explored the role that their tentacles and oral arms may have on their potential swimming performance. We perform comparative in silico experiments to study how tentacle/oral arm number, length, placement, and density affect forward swimming speeds, cost of transport, and fluid mixing. An open source implementation of the immersed boundary method was used (IB2d) to solve the fully coupled fluid–structure interaction problem of an idealized flexible jellyfish bell with poroelastic tentacles/oral arms in a viscous, incompressible fluid. Overall tentacles/oral arms inhibit forward swimming speeds, by appearing to suppress vortex formation. Nonlinear relationships between length and fluid scale (Reynolds Number) as well as tentacle/oral arm number, density, and placement are observed, illustrating that small changes in morphology could result in significant decreases in swimming speeds, in some cases by upwards of 80–90% between cases with or without tentacles/oral arms.